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Abstract 

We investigate a theoretical model applicable to molecular conductors, such 
as TTP and M(dmit)2 salts [M=Ni, Pd, Pt], whose charge carriers orig inate 
from two kinds of molecular orbitals. The model Hamiltonian consists of two 
independent Hubbard chains. The degree of charge transfer between the two 
chains is calculated. The results qualitatively agree with some experimental 
data. 
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I. INTRODUCTION 



In most of the molecular conductors which exhibit metallic properties at low tempera- 
tures, charge carriers originate from molecular orbitals of one kind such as the highest oc- 
cupied molecular orbital (HOMO) of a donor molecule or the lowest unoccupied molecular 
orbital (LUMO) of an acceptor molecule. [1] For example, the charge carriers of /?-(BEDT- 
TTF) 2 I 3 originate from the HOMO's of BEDT-TTF molecules while anions merely work 
as reservoirs of negative charges. In this type of molecular conductors, the chemical compo- 
sition in principle determines the filling factor of the conduction band. 

On the other hand, charge carriers in some molecular conductors originate from two kinds 
of molecular orbitals ( "two-orbital" molecular conductors) . Here we call one of the two kinds 
of molecular orbitals, having the larger orbital energy, as the upper orbital (UO), and the 
other, having the lower orbital energy, as the lower orbital (LO). The energy separation 
between the UO and LO is denoted by A. 

A molecular conductor composed of segregated columns of donor molecules and columns 
of acceptor molecules, such as TTF«TCNQ, may be considered as a two-component type of 
the "two-orbital" molecular conductor. The UO and LO in this type are the LUMO of an 
acceptor molecule and the HOMO of a donor molecule, respectively. 

Another type of the "two-orbital" molecular conductor is a single-component one, where 
the UO and LO come from molecules of one kind. The salt of TTF[Ni(dmit) 2 ]2 is classified 
into this type. The UO and LO in this salt are the HOMO and LUMO of Ni(dmit) 2 , 
respectively. [2] For this type, one can estimate A approximately from the intra- molecular 
optical absorption for the molecule which possesses the UO and LO. The values of A, thus 
determi ned for typical organic donors and acceptors, are listed in Table 1. In the "two- 
orbital" molecular conductors, the charge transfer between the two molecular orbitals is a 
controllable parameter, which one can expect to enrich their properties. In order to reveal 
them, we have studied a two-chain model. 

II. MODEL HAMILTON! AN 

The Hamiltonian which describes two-chained molecular orbitals interacting through the 
on-site Coulomb repulsion is written as 

H = H u + Hi + T + Q, 

H U = AJ2 4a a iv - tun ( a L a i+l* + C "C.) + U uu ^W^W' ( 2A ^ 
i,cr i,cr % 

Hi = -tu ( b t b i+la + CC.) + u u J2 & !tM*W' 

i,cr i 

T = -t ui J2 (4A+1<t + CC.) - U U J2 (btdi+la + CC.) , 
i,cr i,cr 

i,a,<r' 

where a\ a and a ia are the creation and annihilation operator of an electron at the i-th site 
with spin a in the UO; b\ a and b ia are those in the LO; t uu , tu, t u i, U u are transfer integrals; 
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and U uu , Uu, U u \ are the on-site Coulomb repulsions. The model Hamiltonian (1) consists of 
two Hubbard chains (H uu , H u ) and their couplings (T, Q). If the couplings are neglected, 
the problem is reduced to the two Hubbard models which can be solved by Bethe ansatz. 
this approximation may be partially justified by the following considerations: 

First, T is anticipated to be unimportant until the charge transfer from the LO to the 
UO occurs, since t u i, t\ u are smaller than t uu , tu, and negligible compared to A. 

Secondly, 

U 2 

Q = ~y 12 ( a h a n + a li a n + fe W + b li b n) ( 2 - 2 ) 

r 

"2 



12 ( a t a i° + b l b i°) ~ u m 12 ( a \^ a \i a n + b h b iAi b n 1 ■ 



1, a 



where the first term is proportional to the fluctuation of electron density at each site; the 
second term is proportional to the total electron number; and the third term is on-site 
Coulomb repulsion, which simply reduces the effective value of U uu in H u and that of Uu 
in Hi. Consequently, Q can be neglected until the fluctuation of electron density becomes 
dominant. 

Since H u and Hi are now independent, the ground state is obtained by imposing the 
same chemical potential to the two chains. More explicitly one can obtain the ground state 
energy by minimizing E = ^[(UO)*"] + ^[(LO)^] under the condition of N u + N t = N t , 
where ^[(UO)^"] denotes the ground state energy for H u when N u electrons are in the UO, 
and E[(LO) Nl ] denotes the corresponding energy for Hi when iVj electrons are in the LO. 
The total electron number in the system is N t . 

In the following calculation, we use the ground state energy per one site for the Hubbard 
chain. This energy, SHub.( u i n ) = E/N s , was originally derived by Lieb and Wu: [11,12,13] 

rQ 

£Hub(u,n) = —2t / dk cos A; p(k), 
J-Q 

Q 

dk p(k) = N/Ns, (2.3) 

Q 

Q „, ,s 87r „ /4(sin k - smk')\ 



2np(k) = 1 + cos k ( dk' p(k')—R 

J-Q U 
1 f°° i 
R ^ = l + (x + t) 2 



u 

sech(7rt/2) 



-dt, 

where N s , N (< N s ) are the total number of sites and electrons; n = N/N s ; and u = U/t, 
respectively. By solving these equations we can study the range n < 1. For the range n > 1, 
£Hub.( u i n ) i s given by an identical equation, 

£Hub.(u,n) = e Hu b{u, 2-n) + U(n-l), (2.4) 



coupled with Eq. (|2.3| ). Since a closed form of ( |2.3|) is not known, we have obtained following 



functions, Ei(u,n) and 62(11,71) which approximate n): 

£Hub.{u,n) = ei(u,n) (0 < n< 1) 
= e 2 (u,n) (l<n<2) 
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within the ansatz 



£l (u, n) = -A(u) t sm(B(u)n + C(u)n 2 + D(u)n 3 ) 
Ei(u, n) — Ei(u, 2 — n) + U(n — 1); 



where 



A(u 

A(u 

B(u 
C{u 
C{u 
D(u 



2 

- (l + (1 + P±u + P 5 w 2 ) exp(-P 6 w)) (0 < u < 10), 

- + P 13 u~ Pl4 (10<m< 250); 
2/A(u); 

P 7 u/(1 + P 8 u 2 + P 9 u 3 ) + P 10 u (0 < u < 10), 
Pi5U~ Pl6 + P 17 u- Pls (10 < m < 250); 
(tt/2 - B(u)n max - C(u)n 



)/n 



3 . 

max' 



c(«) 



- (l + (1 + P lM + P 2 m 2 ) exp(-P 3 M)) (0 < « < 10), 



(2.5) 



n max (u) = -(1 + P u u~ Pl2 ) (10<u< 250); 

p 1 = -0.08705; P 2 = 0.01366; P 3 = 0.23603; 

P 4 = -0.10512; P 5 = 0.0207; P 6 = 0.28767; 

P 7 = -0.14159; P 8 = 0.06386; P 9 = 0.00582; 

p 10 = -0.00903; Pn = 0.76063; P 12 = 1.03065; 

P 13 = 0.80627; P M = 1.03467; P 15 = -1.3295; 

p 16 = 0.75014; P 17 = 0.12461; P 18 = 0.39471. 

The deviations of ei(u,n) and e 2 ('U,''T') from eHub.( u i n ) are about O.OOlt in the range of 
< u < 250. 



III. RESULTS 

We take U uu = Uu = U and t uu —tu — t for simplicity. Generalization to the other cases 
is straightforward. 

a) (uof-^^LOf-^ 

This is the case corresponding to the charge-transfer salts of TTP 0,5+ .[9] The parameter 
x represents the degree of charge transfer from the LO to the UO. The ground-state energy 
per one site, e(= E/N s ), is expressed by 

£ = (e 2 (U/t, 1.5 + x) + (1.5 + x)A) + e 2 (U/t, 2-x), (3.1) 
= l.W + (1.5 + x)A + s^U/t, 0.5 -x)+ e^U/t, x). 

In what follows, x is determined as the value which minimizes e for a given set oiU/t and 
A/t: 
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> 0, we have, x 



o 

ii) Next, if g 



0. 



x=+0 

de_ 

dx 



< 0, we obtain x as the solution of 



A - 



dei(U/t,n) 



dn 



([//*, n) 



n=0.5- 



dn 



0. 



(3.2) 



Figure 1 shows the degree of charge transfer thus obtained. In the case of TTP salts, [9] 
the value of A/t is considered to be always larger than 2. [10,14] Therefore we conclude that 
the charge transfer from the LO to the UO does not occur in these salts. In fact, properties 
of these salts are consistent with the expectation given by the extended Hiikel calculation 
where only the HOMO is taken into account. [14] 

b) (XJOf +x)Ns (LOf- x)N ° 

This is the case corresponding to the charge-transfer salts of TTP + [9] and [M(dmit)2]~ 
(M=Ni, Pd, Pt). The ground-state energy per site is expressed by 



e = U+(l + x)A + e^U/t, 1 - x) + e^U/t, x). 
The degree of charge transfer, x, is obtained as follows: 



(3.3) 



> 0, we have x = 0. 



ii) Next if g 



x=+Q 



< 0, we obtain x as the solution of 



A 



de x {U/t,n) 



dn 



de^U/t.n) 



n=l— x 



dn 



(3.4) 



Figure 2 shows the degree of charge transfer thus obtained. The ground state of the 
one-dimensional Hubbard model is insulating at half filling. Therefore, the two-chained 
Hubbard model (|2.1|) predicts an insulating state at x — 0, metallic state for x ^ 0, and 
metal-insulator transition on the b oundary of these two regions. 

In the case of TTP + [9] salts, the value of U/t is about 3. [10] Thus , from Fig 2a, one 
may expect that metallic state in these salts if A/t < 3. This condition is not unreasonable. 
In fact, Mori et al recently found the metallic behavior above 160 K in TTM-TTP«I 3 , which 
is a charge-transfer salt of TTM-TTP + . Although they did not consider the possibility of 
the charge transfer between the LO and UO[15], our calculation suggests that such charge 
transfer causes the metallic behavior of this salt above 160 K. 

c) ( UO )(°- 5 +^(LO) (2 ^ 

This is the case corresponding to the charge-transfer salts of [M(dmit)2]°' 5 ~ (M=Ni, Pd, 
Pt). The ground-state energy per one site, e, is expressed as follows: 

e = U(l-x) + (0.5 + x)A + e 1 (U/t,0.5 + x)+e 1 (U/t,x) (0 < x < 0.5), (3.5) 
e = Q.5U+(0.5 + x)A + e 1 (U/t,1.5-x)+e 1 (U/t,x) (0.5 < x < 0.75). (3.6) 

The degree of charge transfer, x, is determined as follows: 
> 0, we have x = 0. 



i) Hi 



x=+0 



5 



< Of* 

x=+0 9x 



=0.5-0' 



we obtain x as the solution of 



U + A + 



de 1 (U/t,n) 



dn 



de^U/t.n) 



n=0.5+x 



dn 



= 0. 



(3.7) 



iv) If 



2=0.5-0 

de 



<o<i 



dx 



a;=0.5+0 



, we have x = 0.5. 

2=0.5+0 

< 0, we obtain x as the solution of 



A - 



de^U/t.n) 



dn 



de^U/t.n) 



n=1.5— x 



dn 



0. 



(3.8) 



Figure 3 shows the degree of charge transfer thus obtained. It should be noted that the 
state x = 0.5 is stable in this model, as can be seen from Fig. 3a. This comes from the fact 
that deHub. {u, n) /dn is discontinuous at n — 1. 

The present model Hamiltonian does not include any contribution from the electron- 
phonon interaction. The interaction is expected to make the state, x = 0.5, more stable by 
forming a gap in the middle of the half-filled UO band (Peierls transition). In this sense, 
the dimerized structure frequently found in Pd(dmit)2 and Pt(dmit)2 salts [16] may be con- 
sidered as a result of Peierls transition. 



d) (V0) xNa (L0f- x)N3 

This is the case corresponding to the "two-orbital" molecular conductors of the two- 
component type, such as TTF«TCNQ. The ground-state energy per one site, e, is expressed 
by 



£ = u(l - x) + xA + 2ei{U/t, x). 
The degree of charge transfer, x, is determined as follows: 



(3.9) 



i) If 



de 

dx 



> 0, we have x 

x=+0 

< < £ 



2 = 1-0 



0. 



, we obtain x as the solution of 



U + A + 



de^U/t.n) 



=1-0 



dn 

< 0, we have x — 1. 



de^U/t.n) 



n=0.5+2 



dn 



= 0. 



(3.10) 



Figure 4 shows the degree of charge transfer thus obtained. This figure exhibits an 
insulating state at x = and x — 1; metallic state for < x < 1. Although a "two- 
orbital" molecular conductor of a single-component type which exhibits metallic behavior 
has not been found, the present model predicts that synthesis of such compound is possible 
by choosing appropriate values of U/t and A/t. 



IV. CONCLUSION 

We found that the on-site coulomb interaction plays quite an important role in the "two- 
orbital" molecular conductors. The effect of this interaction grows the electron density in 
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the UO becomes smaller. The model gives a plausible explanation for the metallic behavior 
of TTM-TTP • I 3 and the strongly dimerized structure frequently observed in Pd(dmit)2 
and Pt(dmit) 2 salts. In spite of its simplicity, this model is versatile and applicable to many 
molecular conductors. 
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TABLES 

s.: determined from an absorption spectrum of solution 

c.: determined from an absorption spectrum of a single crystal 





A(eV) 


UO 


LO 


Ref. 


acceptor 










TCNQ- 


1.3 


LUMO 


HOMO 


[3]'- 


Ni(dmit) 2 


1.1 


LUMO 


HOMO 


[4]-- 


Pd(dmit)^ 


0.9 


LUMO 


HOMO 


[4]-- 


donor 










TTF+ 


2 


HOMO 


2nd HOMO 


[5]-- 


BPDT-TTF+ 


1.7 


HOMO 


2nd HOMO 


[6]-- 


BMDT-TTF+ 


1.4 


HOMO 


2nd HOMO 


[7] c - 


BEDT-TTF+ 


1.4 


HOMO 


2nd HOMO 


[8] c - 


TTP+ [9] 


0.9 


HOMO 


2nd HOMO 


[io] c - 



TABLE I. The energy separation between the UO and LO. The molecular orbitals (LUMO, 
HOMO, and 2nd HOMO) are defined in the neutral state of an isolated molecule. 
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FIGURES 



Fig. 1 (a) The degree of charge transfer (x) obtained for the two-chained orbitals, 
( UO )(i.5+*)jv.( LO )(2-*)jv._ ^ The regions fori = and for o < x < 0.25. 

Fig. 2 (a) The degree of charge transfer (x) obtained for the two-chained orbitals, 
(V0) il+x)N3 (L0) (2 ~ x)Ns . (b)The regions for x = and for < x < 0.5. 

Fig. 3 (a) The degree of charge transfer (x) obtained for the two-chained orbitals, 
( UO )(o.5+*)jv.( LO )(2-*)jv._ ( b ) The regions fo r x = 0; < a; < 0.5; x = 0.5; and 

0.5 < x < 0.75. 

Fig. 4 (a) The degree of charge transfer (x) obtained for the two-chained orbitals, 
(UO) xNs (LO) (2 ~ x)Ns . (b)The regions for x = 0; < x < 1; and x = I. 
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Fig. 1 Tajima et al 
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Fig. 2 Tajima et al 
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Fig. 3 Tajima et al 
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Fig. 4 Tajima et al 



